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ABSTRACT 

The Lick AGN Monitoring Project 201 1 observing campaign was carried out over the course of 1 1 weeks in 
Spring 201 1 . Here we present the first results from this program, a measurement of the broad-line reverberation 
lag in the Seyfert 1 galaxy Mrk 50. Combining our data with supplemental observations obtained prior to the 
start of the main observing campaign, our dataset covers a total duration of 4.5 months. During this time, Mrk 
50 was highly variable, exhibiting a maximum variability amplitude of a factor of ^ 4 in the U -band continuum 
and a factor of 2 in the H/3 line. Using standard cross-correlation techniques, we find that H/? and H7 lag 
the V-band continuum by T^en = 10.64^q|3 and 8.43^[ 2° days, respectively, while the lag of He II A4686 is 
unresolved. The H/3 line exhibits a symmetric velocity-resolved reverberation signature with shorter lags in the 
high-velocity wings than in the line core, consistent with an origin in a broad-line region dominated by orbital 
motion rather than infall or outflow. Assuming a virial normalization factor of / = 5.25, the virial estimate of 
the black hole mass is (3.2 ± 0.5) x 10^ Mq. These observations demonstrate that Mrk 50 is among the most 
promising nearby active galaxies for detailed investigations of broad-line region structure and dynamics. 

Subject headings: galaxies: active — galaxies: individual (Mrk 50) — galaxies: nuclei 
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1. INTRODUCTION 

Observations of broad-line variability in nearby Seyfert 
galaxies play a central role in the interpretation of the demo- 
graphics and cosmological evolution of supermassive black 
holes in active galactic nuclei (AGNs). By measuring the 
time delay between AGN continuum variations and the subse- 
quent response of the broad-line region (BLR) gas, the light- 
travel time across the BLR, and hence the mean BLR ra- 
dius (rBLR), can be directly measured. These reverberation- 
mapping measurements ha ve been carried out for a few dozen 
low-redshift AG Ns (e.g ., iKaspi et all l2000t iPeterson et al] 
l2004t iBentz et al.ll2009al) . and the observed BLR sizes mea- 
sured via H/3 reverberation range from typically a few light- 
days up to several light-months. The BLR size is observed to 
scale with A GN continuum luminosity roughly as tblr oc 
(iBentz et al.t2009b) . and this relationship makes it possible to 
estimate BLR radii from a single observation of an AGN. 

With a direct measurement or estimate of tblr, and as- 
suming virial motion of BLR clouds, it becomes possible to 
estimate the mass of the black hole in an AGN as Mbh = 
/rBLR(Ay)^/G, where AV is the width of the broad line, 
and f is a dimensionles s scaling factor (e .g., lUMch et"an 
[T984t IKaspi etal.1 12OOOI: lOnken et al.1 1200I . This method 
has been used to estimate black hole masses in large sam- 
ples of AGNs out to the highest observed redshifts (for a 
review, see Vestergaard 2011). Currently, nearly all obser- 
vational constraints on the cosmological growth history of 
supermassive black holes depend on masses derived from 
this virial equation. While the assumption of virial mo- 
tion in the BLR has gained support from a variety of con- 
sistency checks (e.g., Peterson & Wandel 2000; OnkeiTetal] 
120041: iNelson et al. 2004), it remains extremely difficult to ob- 
tain direct constraints on the structure and dynamical state of 
the BLR in any individual AGN, and in the absence of such 
constraints, the inferred black hole masses remain subject to 
substantial systematic uncertainty (Krolik 2001). The most 
promising method to examine the kinematics of BLR gas is 
velocity-resolved reverberation mapping, in which the time- 
delay response of emission-line variability relative to contin- 
uum fluctuations can be measured as a function of line-of- 
sight velocity. Velocity-resolved reverberation data can en- 
code a wealth of information about BLR structure on spatial 
scales that are orders of magnitude too srnall to be resolved by 
any other technique (e.g.. IWelsh & Horneill99lt iHorne et all 
12004) . 

Recently, high-cadence observing campaigns have pro- 
duced major improvements in the quality of veloc ity-resolved 
reverb e ration data for the broad B aimer lines (I Bentz et aLl 
I2009at iDennev et aTl 120101) . For the most highly variable 
objects, it is possible to examine the shape of the two- 
dimensional transfer function, which describes the distribu- 
tion ofbroad-line lag response time as a function of veloc- 
ity (IBentz e t al. 2010), and to apply new modeling techniques 
that can directly constrain the BLR geometry and potentially 
test the c ritical assu mption of virial motion (Pancoast et al. 
I2OI lat iBrewer et al.ll201 lb . In order to increase the number 
of objects with data suitable for such analysis, we carried out 
a new reverberation-mapping campaign in Spring 20 11. 

In this Letter, we present our first results for Mrk 50, 
a Seyfert 1 galaxy at redshift z = 0.023. Past observa- 
tion s have found dr amatic variations in its nuclear luminos- 
ity dPastoriza et al.l [T99 1 ), but it has not previously been a 
reverberation-mapping target. During our monitoring pro- 



gram, Mrk 50 exhibited strong variability, and we detect a 
robust reverberation signal and a significant velocity-resolved 
reverberation response across the broad H/3 emission line. 

2. OBSERVATIONS AND REDUCTIONS 
2.1. Photometry 

From 201 1 January 21 until June 13 (all dates are UT), we 
obtained queue-scheduled V-band images using the 0.76 m 
Katzman Autom atic Imagin g Telescope at Lick Observa- 
tory (iFilippenk o et al.l 120011) . the 0.9 m telescope at the 
Brigham Young West Mountain Observatory (WMO), the 
Super-LOTIS 0.6 m telescope at Kitt Peak, the Faulkes Tele- 
scope South at Siding Spring Observatory, and the Palomar 
1.5 m telescope (Cenkoetal. 2006). Exposure times were 
typically 180-300 s. We attempted to observe Mrk 50 on a 
nightly basis, but poor weather and telescope scheduling is- 
sues left some gaps in temporal coverage. 

All images were bias-subtracted and flattened, and cosmic- 
ray hits were removed using the LA-COSMIC routine 
dvan DokkumI [200 1). In order to remove the host galaxy 
and obtain a light curve of the variable AGN flux, we em- 
ployed i mag e subtraction using the HOTPANTS package by 
A. Be c keiPl which is based on the algorithm described by 
lAlardI (120001) . For each telescope, a high-quality template 
image was chosen, and the template was then aligned with 
each night's image and convolved with a spatially varying 
kernel to match the point-spread function of that image. Af- 
ter subtracting the scaled template image, the variable AGN 
flux is left as a point source in the subtracted image, allow- 
ing for aperture photometry using the IRAF DAOPHOT pack- 
age. The photometric aperture radius used for each telescope 
was set to match the average point-source full width at half- 
maximum intensity (FWHM) for images from that telescope. 
Light curves were constructed separately for each telescope, 
and each was then normalized to the WMO light curve by de- 
termining an average flux scaling factor based on nights when 
Mrk 50 was observed at both locations. We find that the im- 
age subtraction yields a better-quality light curve than simple 
aperture photometry, and provides closer agreement between 
the light curves obtained with different telescopes. 

The overall flux scale was calibrated using observations of 
Landolt (1992) standard stars observed during a few photo- 
metric nights at WMO. Observations taken within 6 hr of one 
another were combined by taking a weighted average. The 
final V fight curve is illustrated in Figure [T] 

2.2. Spectroscopy 

Our campaign at the Lick Shane 3 m telescope consisted 
of 69 nights allocated between 2011 March 27 and June 13, 
during which time we obs erved Mrk 50 on 41 n ights using 
the Kast dual spectrograph dMiller & StonelfT993l) . In this pa- 
per, we discuss only measurements from the blue arm of the 
spectrograph, where we used a 600 lines mm"' grism over 

3440-5520 A at a scale of 1.0 A pixel"'. All observations 
were done with a 4"-wide slit oriented at PA = 180°. Standard 
calibration frames including arc lamps and dome flats were 
observed each afternoon, and flux standards were observed 
during twilight. The exposure time for Mrk 50 was normally 
2 X 1200 s. Additionally, we observed Mrk 50 on 18 nights 
that were allocated to other observing programs, beginning 
on January 26. All observations used the same 600-line grism 

http://www.astro. Washington. edu/users/becker/c_softwai'e. html 
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Fig. 1 . — Mrk 50 light cui-ves for the V band, the J/,-band continuum mea- 
sured from the spectra over 3550-3800 A, and the H/3, H7, and He II emis- 
sion lines. The V-band panel shows a difference imaging light curve illus- 
trating changes in flux relative to the mean. Units for the V and Us bands are 
10"'^ erg cm"^ s"' A"', and units for broad-line fluxes are 10"'^ erg cm"^ 
s"' . The arrow in the H/3 light curve marks the start of the dedicated Lick 
observing campaign. 



and 4" slit, but the wavelength coverage was slightly different 
each time. The exposure for these observations was typically 
900 s. 

Spectroscopic reductions and calibrations followed stan- 
dard methods implemented in IRAF and IDL. A large extrac- 
tion width of 10"3 was used in order to accommodate the full 
extent of the AGN spatial profiles observed on nights with 
very poor seeing. Error spectra were extracted and propa- 
gated through the full sequence of calibrations, and for each 
night the two exposures of Mrk 50 were averaged together. In 
the region 4600^800 A, the median signal-to-noise ratio per 
pixel is 75. 

3. SPECTROSCOPIC DATA ANALYSIS 

The reduced spectra were first normalized to a 
uniform flux scale by emplo ying the procedure of 
Ivan Groningen & WandersI (119921) . This method applies 
a flux scaling factor, a linear wavelength shift, and a Gaussian 
convolution to each spectrum in order to minimize the 
residuals between the data and a reference spectrum con- 
structed from several of the best-quality nights. The scaling 
is determined using a wavelength range containing the [O III] 
A5007 line, which is assumed to have constant flux. 

Since the spectra contain a substantial amount of host- 
galaxy starlight, we implemented a simple continuum sub- 
traction routine to produce a cleaner measurement of the 
broad-line profiles an d fluxe s, similar to the procedure de- 
scribed by iPark et al.l (1201 II) . Over the wavelength range 
4300-5400 A, each scaled spectrum was fitted with a model 
consisting of a power-law featureless continuum, an Fe 11 tem- 




4600 4800 5000 5200 5400 
Observed Wavelength (A) 

Fig. 2. — Mean and rms spectra. In each panel, the upper spectrum is 
constructed from the set of scaled spectra of Mrk 50, and the lower spectrum 
is constructed from the continuum-subtracted, scaled spectra. 



plate from iBoroson & GreenI (119921) broadened by convolu- 
tion with a Gaussian kernel, several emission-line compo- 
nents r epresented by either Gaussian s or Gauss-Hermite func- 
tions ( van der Marel & Franxl 11993). and a starlight model 
consisting of simple stellar pop ulation models at solar metal- 



jpul; 

licity from ' Bruzual & CharlotI (l2003l) which were broadened 
by convolution with a Gaussian kernel. Additionally, a 
foreground extinction was applied to the model spectrum. 
The model fit was optimized by minimization using a 
Levenberg-Marquardt technique (Markwardt 2009). Then, 
the best-fitting model components representing the starlight 
and nonstellar continuum were subtracted from the data, leav- 
ing a pure emission-line spectrum. For the starlight model, 
we obtained a good fit using an 11 Gyr-old population as the 
dominant component, and adding contributions from younger 
populations did not significantly improve the fit. The median 
starlight fraction at Arest = 5 100 A is 41 %, and Fe II contributes 
just ^ 2-3% of the continuum flux density at 4600 A. 

In prior reverberation work, starlight subtraction has not 
generally been applied, and the broad-line light curves have 
typically been measured by subtracting a local, linear contin- 
uum fitted to wavelength regions on either side of an emission 
line. We find that our continuum subtraction procedure pro- 
vides a more accurate removal of the continuum shape under- 
lying each emission line and a much better light curve for the 
weak He II A4686 line. Figure |2] shows two versions of the 
mean and root-mean-square (rms) spectra: one constructed 
from the set of scaled spectra, and one constructed from the 
set of continuum-subtracted, scaled spectra. This illustrates 
the utility of the continuum subtraction procedure for remov- 
ing the stellar absorption features and continuum undulations. 

The broad-line light curves were then measured by integrat- 
ing the flux over fixed wavelength ranges in the continuum- 
subtracted spectra: 4890-5050 A for H/3, 4370-4510 A for 
H7 (also including the narrow [O III] A4363 line), and 4730- 
4870 A for He II. Emission-line light curves are shown in Fig- 
ure [U along with the continuum flux density measured from 
the scaled spectra over 3550-3800 A (which we refer to as the 
Us band). The continuum and broad lines were highly variable 
during the monitoring period: over a 30-day span the f/j-band 
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Fig. 3. — Cross-correlation functions for H/9, H7, and He 11 against the 
AGN continuum flux, the autocorrelation function of the V-band continuum, 
and the cross-correlation of the continuum bands Us vs. V . 

continuum dropped by a factor of ^ 4, and the H/3 line re- 
sponded with a factor of ~ 2 decline. For observations taken 
before the start of our dedicated campaign on March 27, the 
scatter in the spectroscopic light curves is relatively high. We 
attribute this to the nightly differences in observational setup 
and calibrations used by each observing team, and to the very 
poor observing conditions during some winter nights. 

To estimate black hole masses from reverberation data, the 
second moment of the H/3 line in the rms s pectrum is mos t 
often used as the measure of line width (Pete rson et al.ll2004t) . 
We find a{\\fi^^^ = 1740 ± 101 km s"', where the measure- 
ment uncertainty is determined th rough a bootstrap resam- 
pling procedure dBentz et"an i2009a). and the line width has 
been corrected for instrumental broadening of ainst ~ 133 
km s"' following lBarth et alj (120111) . 

4. REVERBERATION LAG MEASUREMENTS 

In order to measure the cross-correlation function (CCF) 
for unevenly sampled time series, we employ the interpola- 
tion cross-correlation function methodolo gy and Monte Carlo 
error a nalysis techniques described by 'Gaskell & Peterson' 
(ll987h . lWhite & PetersonI (Il994t) . and|Peterson et al. (2004); 
these methods have been employed in the majority of re- 
cent reverberation- mapping work (e.g., Bentz et al. 2009a; 
iDennev et al.ll2010l) . We measured the cross-correlations of 
H/3, H7, and He II against the V light curve, and also of H/3 
against the Us light curve. The CCFs were computed from 
-20 to +40 days in increments of 0.25 days. Table[T]lists two 
measures of the lag: Tpeak, which is the lag at the peak of the 
CCF, and Teen, the centroid of the CCF for all points above 
80% of the peak value (P eterson et al.,.2004i) . Figure |3] illus- 
trates the CCF measurements. 

The H/3 and H7 lines have similar lag times of Teen = 10.64 
and 8.43 days, respectively, but for He II we find that both Teen 
and Tpeak are consistent with zero lag, indicating a very com- 
pact size for the inner, high-ionization portion of the BLR. 
The faster response time for He II is apparent in the light 
curves, particularly in the steep drop beginning in mid-April 
(HJD w 2455670). The Balmer lines, in contrast, show a 
much more gradual and delayed decline in response to the 
falling continuum flux. The flattened peak and asymmetry 
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Note. — All lags are given in the observed frame. 

of the H/3 CCF in comparison with the continuum autocorre- 
lation function (ACF) suggests that the Balmer-line emitting 
zone of the BLR covers a fairly large radial extent, and this 
interpretation is supported by the velocity-resolved measure- 
ments described below. We also measured the CCF between 
the Us and V bands in order to search for evidence of rever- 
beration due to Balmer continuum emission from the BLR 
dMaoz et al. 1993; Korista & Goad 2001), but no significant 
lag was found. 

Our Mrk 50 dataset, which benefits from high-amplitude 
variability and high-cadence sampling, presents an excellent 
case study for velocity-resolved variability. Light curves were 
measured for seven separate velocity segments across the 
width of the H/3 line, and each segment light curve was cross- 
correlated against the V light curve. Figure |4] illustrates the 
lag (Teen) as a function of velocity across the H/3 line, reveal- 
ing a roughly symmetric trend of shorter lags in the line wings 
and longer lags in the core, with a ~ 10 day difference in re- 
sponse time between the core and wings. This pattern, which 
resembles the symmetric H/3 lag res ponse seen in some other 
AGNs inclu ding Mrk 1 10 jKoUatschnv & Bischoff 2002) and 
NGC 5548 (IDennev et al J 1201 01) . is qualitatively consistent 
with predictions for BLR clouds in cir cular orbits in th e Kep- 
lerian potential of the black hole (e.g.. lWelsh & Horn3 [l991), 
with the highest line-of-sight velocities originating from gas 
located close to the black hole. In contrast, a BLR dominated 
by either radial infall or outflow would result in an asymmet- 
ric velocity-lag map, with inflow producing longer lags on the 
blueshifted side of the line, and outflow producing longer lags 
on the red side. Suc h signatures of radial motion have bee n 
seen in a few objects (iBentz et alJ2009al ; lDennev et al.l l2010l). 
but the current sample of AGNs with velocity-resolved data 
of this quality is still too small to examine the distribution of 
different BLR kinematic states. 

5. BLACK HOLE MASS ESTIMATE 

Following iPeterson et al.l ( 120041) . we compute the H/3 
"virial product" [defined as VP = rBLR(Ay)^/G] by using 
Ay = (7(H/3rms) and tblr = cTeen, whcrc Teen has been cor- 
rected to the AGN rest frame. For Mrk 50, tblr = 10.40;o 
It-days, and VP = (6.2 ±0.9) x 10^ Mq. 

While the virial product is a straightforward combination 
of measured quantities, the relationship between virial prod- 
uct and black hole mass is more indirect and uncertain. Most 
recent reverberation work has adopted a normalization of 
the virial mass scale based on the assumption that AGNs as 
a whole fall on the same Mbh - <^ relation as nearby inac- 
tive galaxies dOnken et al. 2004; W oo et aiJ l2010). Using the 
A^BH - c relation derived by Giiltekin et aljl2009 ) as the lo- 
cal reference, this implies a mean value of log/ = 0.72;!;Q 'j'Q, 
and the reverberation masses determined in this way show 
an intrinsic scatter of 0.44 dex about the Mbh - f relation 
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Fig. 4. — Velocity-resolved reverberation in the H/3 line. The upper panel 
shows the mean lag measured for each velocity segment of the broad H/3 
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(I Woo et al.ll20Tol) . Adopting the Woo et al. value of / = 5.25 
implies that Mbh = (3 .2 ± 0. 5)x lO^M^ forMrk 50. For con- 
sistency wit h recent work (iPeterson et alJ l2004t iBentz et aP 
l2009atlDen ney et al. 2010!), the quoted uncertainty includes 
only the statistical error on the virial product, but not the un- 
certainty resulting from the choice of a specific / value; the 
true uncertainty in Mbh is dominated by the uncertainty in 
/. Masses derived in this way depend on the assumption of 
virial motion, the assumption that AGNs should fall on the 
same Mbh ~ relation as quiescent galaxies, and the adop- 
tion of a specific Mbh-o" relation as the local baseline. Dif- 
fering assessments of the form of the local Mbh-ct relation, 
particularly at low masses, can potentially affect the normal- 
ization of the AGN virial mass scale at th e factor of 2 level 
(iGreeneet al.ll2010l: iGraham et alll2011h . Additionally, there 
is ongoing debate over the possible effect of radiation pres- 
sure on BLR clouds and its impact on the i nferred masses 
dMarconi et al.ll2008l: iNetzer & MarzianillloTQ) . Thus, while 
the mean value of / can currently be determined to ^ 20% 
precision within the context of a specific set of assumptions 
about the Mbh-ct relation (Woo et al. 2010; Graham et al. 
1201 Ih . the actual uncertainty in the overall AGN mass scale 



remains significantly larger and is difficult to quantify, and 
the observed 0.44 dex scatter in the AGN Mbh - c relation 
must reflect, at least in part, the dispersion of the true / 
values for individual AGNs. Resolving these issues will re- 
quire enlarging the number of targets having high-quality 
reverberation-mapping data, and the application of new ap- 
proaches t o extract dynarnical information from the observa- 
tions (e.g., IHncoiiLSiDliniil IlliwEiEiD'^MB 

Mrk 50 has an early-type morphology, and MaUca net alj 
([T998h classify it as an SO galaxy based on Hubble Space 
Telescope imaging. The only published measurement of the 
stellar velocity dispersion of Mrk 50 is (T* = 7 8 ± 15 km s~' 
based on a Sloan Digital Sky Survey spectrum (iGreene & Hd 
120061) . However, a new measurement from Keck LRIS data 
gives 109 ± 14 km s"' (Harris et al., in preparation), and we 
consider this to be more reliable than the SDSS measurement. 
The scaling relations calibrated bv lGiiltekin et al.l (12009 ) then 
imply an expected Mbh ~ (0.6- 1.7) x 10^ Mq from the 
general Mbh - c relation including ellipticals and spirals, or 
(0.9-2.5) X 10^ M0 based on the Mbh -cr relation fitted toel- 
lipticals only. Our reverberation-based mass is slightly higher 
than these values, but Mrk 50 still lies well within the ^ 0.4 
dex scatter of the AGN Mbh - relation ( Woo et al.,.2010i) . 

6. CONCLUSIONS AND FUTURE WORK 

From our Spring 2011 monitoring campaign, we have ob- 
tained very high-quality reverberation mapping of Mrk 50. 
This is one of just a few nearby AGNs in which strong 
velocity-resolved lag signatures have been detected, and it is 
now among the most promising targets for detailed studies of 
BLR geometry and kinematics. 

Our long observing campaign produced a large and very 
rich dataset, and this paper presents just one set of early re- 
sults from this program. Future papers will include detailed 
descriptions of the data-analysis procedures and results for 
the entire observed sample. A major goal of our project is to 
exploit the potential of velocity-resolved reverberation map- 
ping to elucidate the structure of the BLR and to derive black 
hole m asses directly, and an upcoming paper (Pancoa st et alj 
1201 Ib l) will describe new modeling of our Mrk 50 data. The 
black hole mass determined from dynamical modeling is con- 
sistent with the simple virial estimate present ed here, and 
Mrk 50 is now the second object (after Arp 151; lBrewer et alj 
1201 II) to show agreement between the two approaches. 
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